Phonon-assisted relaxation between hole states in quantum dot molecules 
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We study theoretically phonon-assisted relaxation and inelastic tunneling of holes in a double 
quantum dot. We derive hole states and relaxation rates from k p Hamiltonians and show that 
there is a finite distance between the dots where lifetimes of hole states are very long which is related 
to vanishing tunnel coupling. We show also that the light hole admixture to hole states can affect 
the hole relaxation rates even though its magnitude is very small. 
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Systems composed of two stacked quantum dots are in- 
teresting due to their wide spectrum of non-trivial phys- 
ical properties. In closely spaced double quantum dot 
structures, carrier states can be strongly delocalizecP 
(like in a chemical covalent bond in molecules), hence 
these systems are often called quantum dot molecules 
(QDMs). They are promising for many possible applica- 
tions, e.g., in quantum computing and nanoelectronics . 
Holes in QDMs can be particularly interesting due to 
their possible applications in quantum computing. Be- 
cause of the long hole spin lifetime^ and long coherence 
times, holes has been proposed as realization of qubiiP^. 
Hole states can be optically controlled using picoseconds 
optical pulses^. Furthermore, recent results show that 
hole spin states can be prepared with high fidelity^. How- 
ever, phonon-related processes are inevitable in a crystal 
environment and may limit the feasibility of quantum 
control in these systems. 

The properties of the hole states are nontrivial as a 
result of subband mixing. In particular, recent theoreti- 
cal predictions^^ an( j experiments^ indicate that hole 
states in QDMs show an unusual behavior: the ground 
state of the hole becomes antibonding above a certain 
critical distance between the dots. Thus, there is a dis- 
tance where degeneracy between the bonding and anti- 
bonding molecular states occurs. It can be expected that 
this behavior will be reflected in the phonon-assisted tun- 
neling rates between the dots, leading to qualitatively 
new behavior that was not present in the widely stud- 
ied electron case^HUl These effects result from the non- 
trivial structure of hole states (eg. subband mixing) and 
cannot be revealed in simple models of hole states. Hole 
relaxation rates have been recently measured^ but only 
a simple model was used to interpret the results, where 
band mixing effects were not taken account. 

In this work, we present theoretical results for hole 
relaxation rates in a self-assembled QDM. In view of 
the important role of subband mixing for hole states in 
QDMiP^HEfl we propose an approach based on multi- 
band k ■ p method with the hole-phonon coupling de- 
scribed by the Bir-Pikus strain Hamiltonian. Our results 
show that for a certain finite distance (where the tunnel 
splitting vanishes) , phonon-assisted relaxation is dramat- 
ically slowed down. We also prove that even small (below 
2%) light hole admixture to hole states can strongly af- 



fect the deformation potential coupling between carriers 
and phonons (up to 25% of the corresponding relaxation 
rate). 

The paper is organized as follows. First, we define 
the model of the system under consideration. Next, we 
present the method used to calculate hole states. Subse- 
quently, phonon-assisted relaxation of hole states is dis- 
cussed. 

We consider an axially symmetric system formed by 
two geometrically identical self-assembled InAs dots in a 
GaAs matrix. The shape of both the dots is modeled as 
a spherical segment with the base radius of 10 nm and 
the height 3.7 nm. Dots are placed on a wetting layer 
with thickness 0.6 nm. The dots are separated by a dis- 
tance D (base to base) which is a variable parameter of 
our model. A diffusion layer of a very small thickness 
0.3 nm is included at the contact between the two ma- 
terials, where the material composition changes linearly. 
The system is placed in an axial electric field. The strain 
due to lattice mismatch is represented by the strain ten- 
sor e and calculated by minimizing the elastic energy of 
the system in the continuum elasticity approximation^!. 

We use the 8-band Kane Hamiltonian with Bir-Pikus 
strain dependent termsP 3 24 which introduces dependence 
on position via the position dependent strain fields and 
parameter discontinuities at the InAs/ GaAs border. We 
apply Burt-Foreman operator ordering^ 25 ! 26 !. U sing the 
quasi degenerated Lowdin perturbation methocpSEZl we 
reduce the problem to the 4x4 part describing the heavy 
hole (hh) and light hole (lh) subbands of the valence 
band. For all the elements of this part we calculate the 
perturbative corrections resulting from the coupling to 
the conduction band up to order fc 2 . Corrections of order 
k 2 for diagonal elements are expressed in terms of (po- 
sition dependent) effective masses. The corrections from 
the spin-orbit subbands are of the order k 4 and therefore 
are neglected. 

The matrix representation of the effective valence band 
Hamiltonian in the basis {|hh f), |lh f), |lh 1), |hh |)}, 
where j" and J, represent the projection of the total angu- 
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lar momentum (±| for hh and ±i for lh), is 
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where (in cylindrical coordinates) 
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Here e z is an electric field, 7 ^ are modified Luttinger pa- 
rameters^, the operators k± in cylindrical coordinates 
are given by 



fc ± = e ±_ a z a" 



77io is the free electron mass, effective masses are defined 
as 
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where 

E s ,hb(p, z) = E c (p, z) - E hh (p, z), 
Eg,ib(p, z) = E c (p, z) - E lh (p, z). 

Position-dependent band edges E c (p,z), £"hh(p, z) and 
-f?ih(p, z) are given by 

E c (p, z) = E c0 + a c Tr{e}, 
Ehh(p, z) = E v0 - a„Tr{e} - b v [e zz - 0.5(e xx + e yy )], 
z) = E v0 - a v Tr{i} + b v [e zz - 0.5(e xx + e yy )], 



Band structure parameters _E c o 

-Evo 
E P 

Modified Luttinger parameters 7 J 

72 
73 

a c 
a v 
6 
ci 
c t 



0.95 eV 

-0.57 eV -0.42 eV 

25.7 eV 22.2 eV 



Deformation potentials 



1.34 
-0.57 
0.062 
-9.3 eV 
0.7 eV 
-2.0 eV 



1.98 

-0.44 

0.48 

-6.66 eV 
0.66 eV 
-1.8 eV 



Speed of sound - longitudinal 

- transverse 
Crystal density 
Piezoelectric constant 
Relative dielectric constant 



5150 m/s 
2800 m/s 
5300 kg/m 3 
-0.16 C/m 2 
12.9 



TABLE I: Material parameters used in the calculations. 



where a v and b v are the valence band deformation po- 
tentials and a c is the conduction band deformation po- 
tential. The values of the material parameters are listed 
in Table 3PI The band offsets assumed here are close to 
the values in Ref. [25] 

The Hamiltonian describing the hole-phonon coupling 

;J3Q1 
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where tjj n is the four-component eigenfunction the el- 
ements of which are related to the valence subbands, 
Hb-p is the valence part (4x4) of the Bir-Pikus Hamilto- 
nian describing the deformation potential (DP) coupling, 
Hpe — IV (r) is the part of the Hamiltonian describ- 
ing piezoelectric (PE) coupling and / is the 4x4 identity 
matrix. The piezoelectric potential is given by V(r) — 
—i(d- e)||/(eoe r ) where d is the piezoelectric tensor, e is 
the vacuum permittivity and e r is the relative dielectric 
constant of GaAs. In the zincblende structure, the only 
non-zero components of d are d xyz = d yzx = d zxy = d. 
Because of the orthogonality and identical parity of Bloch 
functions for different valence subbands, the lowest order 
inter-subband contribution from the piezoelectric cou- 
pling contains an additional factor on the order of (qa) 2 , 
where a is the lattice constant. This term not only qual- 
itatively affects the small q behavior of the coupling con- 
stant but also is quantitatively small for low values of 
q. Therefore, for long wavelength phonons, the inter- 
subband elements of the piezoelectric potential are much 
smaller (by 2 orders of magnitude in GaAs) than those 
resulting from the DP coupling. We take into account 
only acoustic phonons because the inter-level energy dis- 
tance in our structure is much smaller than the optical 
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phonon energy. We assume linear dispersion of acoustic 
phonons. 

In order to describe the hole states and account 
for the subband mixing effects we propose a multi- 
subband gene ralization of the 'adiabatic' separation of 
variables' 14 * 31 ^. For each subband, we solve numerically 
the one-dimensional equation along the growth (z) direc- 
tion for each value of the radial coordinate p, 
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where a is a subband index. We find the two lowest 
solutions to this equation for the hh and lh subbands 
and obtain Xhh,n{p, z),xih,n{p, z), n — 0, 1 as well as 
the corresponding p-dependent eigenvalues which can be 
interpreted as effective potentials for the radial prob- 
lem. In the next step, we apply the Ritz variational 
method, looking for the stationary points of the func- 
tional F[ip] — (tp\H\i/}). We use the class of normalized 
trial functions 



ip(p, Z, (j>) 
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where M is the projection of the angular momentum on 
the system axis and the four components refer to the 
four valence subbands. Finally, we numerically minimize 
the functional F[ip] on a grid and obtain the two lowest 
eigenfunctions of the system. 

The obtained lowest eigenstates have a small lh admix- 
ture (from 0.5% to 2% contribution to the total state) due 
to a very weak confinement of light holes in our structure. 
In order to characterize the tunnel coupling between the 
hole states in the two dots, we study the anticrossing be- 
tween the two lowest states in the QDM as the electric 
field is scanned through resonance. The effective tunnel 
coupling parameter is defined as half of the anticrossing 
width, with a positive (negative) sign assigned in the case 
of bonding (antibonding) ground state. The tunnel cou- 
pling parameter t is shown in Fig. |TJb) , as a function 
of the distance between the dots D. The transition from 
bonding to antibonding ground state, which occurs about 
D = 10.23 nm in our structure, corresponds to vanishing 
tunnel coupling between the dots. 

We calculate phonon-assisted relaxation rates using 
the Fermi golden rule^l with the Hamiltonian ([2| . The 
total phonon-assisted relaxation rate for two different 
distances between the dots (D = 7 nm and 12 nm) 
are shown in Fig. |2ja,b). The rates are calculated as 
a function of the axial electric field at three tempera- 
tures (T = 0, 20 K and 40 K). The contributions from 
the DP and PE couplings to the total relaxation rate 
are shown in the Fig. |2jc-f). For closely spaced dots 
(D = 7 nm) the energy splitting 2t is high. Then, relax- 
ation rates from DP and PE couplings are comparable, 
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FIG. 1: (Color online) (a) Energy of the two lowest hole states 
for D — 12 nm. The absolute value of the tunnel coupling 
parameter is equal to half of the energy splitting at the res- 
onance, (b) The tunnel coupling parameter as a function of 
the distance between the dots. The point of t = is where 
the bonding-antibonding transition of the ground state takes 
place. 



although relaxation rate from DP slightly dominates. On 
the contrary, for larger distances, the PE interaction is 
much stronger. This behavior is similar to the electron 
case^. The DP coupling dominates for large energy split- 
tings because this coupling is isotropic and involves LA 
phonons which have higher energies, while the PE cou- 
pling is anisotropic and is suppressed for emission along 
the z direction which is preferred at high energies^. On 
the other hand, for small splittings, the PE coupling be- 
comes realatively stronger because of its strong 1/q de- 
pendence at q — » 0. 

In order to study the overall dependence of the relax- 
ation rate on the inter-dot separation, we calculate the 
maximum (over field magnitudes) phonon-assisted relax- 
ation rate as a function of D. The total maximum relax- 
ation rate as well as the contributions from the DP and 
PE couplings are shown in Fig. [3^a-b). For D — > both 
the relaxation rates drop down because the energy split- 
ting is large and the density of phonon states at very high 
frequencies is low. On the other hand, for D — > oo both 
the relaxation rates decay due to vanishing overlap be- 
tween the wavefunctions. One can see that for distances 
near D = 10.23 nm (where t — 0) the maximum phonon- 
assisted relaxation rate drops down to t -1 = 0.0798 ns _1 
and r _1 = 1.04 ns _1 for the PE and DP couplings re- 
spectively, that is, by two orders of magnitude compared 
to the highest values. This results from an extremely 
narrow energy splitting in view of low density of states 
of low frequency phonons. Clearly, the relaxation at the 
critical distance remains very slow at any electric field. 

We have also verified the accuracy of the single band 
approximation. We have made calculation of the phonon 
assisted relaxation rate for D = 10.23 nm (Fig. |4ja,b)) 
and for D = 5 nm (Fig. |4^c,d)) and compared the results 
for the DP (Fig. ga^)) and PE couplings (Fig.gb.d)). 
For the DP coupling at D = 10.23 nm, the difference 
between the full model and the single band approxima- 
tion is considerable. This can be interpreted as follows. 
Near the point t = 0, the two lowest eigenfunctions are 
strongly localized on different dots and the overlap of 
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FIG. 4: (Color online) (a,c) Contribution to the relaxation 
rate due to the DP coupling. Results from 4-band model 
(red solid line) and from the single band approximation (blue 
dashed line). (b,d) The same comparison for the PE coupling. 



FIG. 2: (Color online) (a-b) Total phonon assisted relaxation 
rate for two different distances between dots at T = K (red 
solid line), 20 K (blue dashed line), and 40 K (green dotted 
line). (c,d) Contributions to the phonon-induced relaxation 
rate due to DP coupling. (e,f) Contributions from PE cou- 
pling. 
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FIG. 3: (Color online) (a) The total maximum relaxation rate 
(blue dashed line) and the contribution from the PE coupling 
(red solid line), (b) Contribution to the total relaxation rate 
from the DP coupling. 



their hh components is very small. However, lh compo- 
nents of these wavefunctions are still delocalized due to 
the shallow confinement in the lh subband. This leads 
to a strong coupling between hh component of the first 
state and the lh component of the second one via the 
inter-subband elements of the strain (Bir-Pikus) Hamil- 
tonian. On the other hand, for PE interaction, the mis- 
match between the models is very small. The reason 
is the absence of inter-subband terms for PE coupling. 
We can also see that for both the DP and PE cases the 
relaxation rate vanishes at a certain electric field corre- 
sponding to the exact resonance. At D = 5 nm (Fig. Eb) 
for DP coupling, the discrepancy between the two mod- 



els is much smaller than for D — 10.23 nm. For stronger 
tunnel coupling (smaller D), the resonance area, where 
the wavefunctions are delocalized, extends over a wide 
range of electric fields. In consequence, the interband 
coupling is not so important in this case. For even larger 
values of D, localization is strong, but lh contribution 
to the states decreases, hence the contribution from the 
interband coupling is also small. 

In summary we have studied phonon-assisted relax- 
ation between the two lowest hole states in a QDM, fo- 
cusing on the role of subband mixing and the bonding- 
antibonding transition of the ground state. We have 
shown that the relaxation rate is very small at the crit- 
ical distance (degeneracy point), where the character of 
the ground state wavefunction changes. Our findings 
are consistent with the general features of experimental 
observations^, where hole lifetimes on the order of a few 
nanoseconds were observed (the geometry of the system 
studied in that work was, however, different from that 
considered here). We have also investigated the influ- 
ence of the lh admixture on the relaxation rate. We have 
shown that near the degeneracy point, the subband mix- 
ing gives an important correction to phonon-induced re- 
laxation via DP channel, up to 25% in some cases, in spite 
of a very small contribution of the lh subband to the hole 
state (below 2%). Our results show that phonon-assisted 
hole relaxation and tunneling in QDMs are with some re- 
spects qualitatively different from the electron case. The 
complete description of the hole-phonon kinetics turns 
out to be impossible without allowing for subband mix- 
ing in contrast to the electron case, where simple models 
are able to essentially correctly account for the relaxation 
rates. 
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